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Abstract 

We present a systematic study of the geometries, energetics, electronic 
structure and bonding in various Al-Li clusters viz. Al n Li n (n=l-ll), ALJ, 
Alg - , Ak^Li, A^Li - , and ALgLig using Born-Oppenheimer molecular dynamics 
within the framework of density functional theory. The growth pattern in 
these clusters is divided in two broad categories: the first one consisting of a 
bent rhombus of AI2LI2 (n=2-4) and the second one consisting of a pentagonal 
ring (n=7-9,ll). A substantial charge transfer is seen in nearly all clusters 
with the exception of Al2Li2 where the charge transfer is not significant. In 
clusters with more than six Al atoms, the eigenvalue spectrum is divided 
into two groups: a lower group of jellium-like states and a higher group of 
localized bonds formed through the interaction of the p electrons on each of 
the Al atoms. Finally, we have discussed the tetravalent behavior of Al atoms 
arising due to a charge transfer from the Li atoms to the Al atoms. 
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I. INTRODUCTION 



The discovery of C 60 buckminsterfullerene [1], observation of magicity in metallic clus- 
ters [2], unusual thermodynamical properties of clusters of atoms like Na [3], Sn [4] etc., 
has given rise to an explosive growth of research in the field of cluster physics. Extensive 
theoretical and experimental [5-7] work has been carried out to understand the physical 
and chemical properties of clusters like size-evolutionary pattern in cluster geometry, ther- 
modynamics, vibrational and rotational properties, optical properties, electronic structure, 
and bonding, chemical reactivity, as well as the emergence of the bulk properties with the 
increase in cluster size. 

Recently observed aromaticity [8] and antiaromaticity [9-12] in inorganic molecules, and 
in particular Al-Li clusters, has stimulated further research in these systems. These prop- 
erties are well-known and important in organic chemistry [13]. However, recently the first 
experimental and theoretical evidence of aromaticity in all-metal system viz. Al^~ and 
MAlJ (M = Li, Na, Cu) was reported by Wang and coworkers [8]. They found the Al| 
dianion to be square-planar possessing two 7r-electrons, thus conforming the structural cri- 
terion for aromaticity. The concept of antiaromaticity, introduced by Breslow [14] et al, 
corresponds to the destabilization seen in the cyclic systems with An 7r-electrons. Although 
such molecules have not yet been observed except in organic chemistry, calculations reported 
by Shetty [9,11] et al. and others [12] indicate the Al 4 Li 4 cluster to be antiaromatic. 

Aluminum-lithium clusters, which we study in this paper, are also interesting due to 
their unusual structural and bonding properties that are quite different from those of the 
pure clusters of the constituent elements. Both pure Al and Li clusters have been extensively 
studied using the spherical jellium model (SJM) [2]. However, it is observed that although 
the SJM has been quite successful in describing the gross electronic structure and the sta- 
bility of alkali metal atom clusters [2,15], the Al clusters present an interesting contrast. 
Rao and Jena [16] have reported a comprehensive study of the aluminum cluster using ab 
initio density functional theory. They found that in small Al clusters, the effective valency 
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of Al is one [16] due to a large energy gap of about 5eV separating the 3s 2 and 3p x orbitals. 
However, Rao [17] et al. as well as Dhavale [18] et al. did not find any signature of such 
monovalent behavior, contrary to the expectations. Investigation by Cheng [19] et al. shows 
that a single Al atom in Li clusters introduces a localized bond between the Al impurity 
and the Li host. They also found a magic cluster AlLi 5 , and suggested that Al n Li 5n can 
aggregate at least for some values of n > 2 and may exhibit properties characteristic to an 
assembly of such clusters from the AlLi 5 subunits. Akola [20] et al, however, found that 
this idea does not apply after n = 2. Further, their investigation on bonding in Li-rich Al-Li 
clusters show a substantial charge transfer from Li atom and nearby Al atom, strengthening 
the ionic Al-Li bond, while the Al-Al bond gained a more covalent nature. Such a charge 
transfer has been observed by Kumar [21] in several mixed Al-Li clusters. He found a layered 
Al 10 Li 8 compound to be magic with electronic and geometric shell closing. This was the 
first instance where a shell closure in s — p bonded metal cluster was found to occur at 38 
valence electrons. Note that the magic numbers for alkali metal atom clusters [2] are 2, 8, 
18, 20, 40,..., whereas, some of the magic Al clusters are Aly , Alf, Al^, Alj~ 3 , etc [16,22]. 

The Al-Li bulk is especially interesting due to the fact that it forms a stable alloy at 
over the wide range of concentration. However, the most stable intermetallic B32 phase 
corresponds to the 50:50 concentration [23]. In this B32 phase, a mixture of covalent and 
ionic bonding is seen due to a remarkable charge transfer from the Li atoms to the Al 
atoms. Moreover, the close resemblance of the density of states (DOS) (without any band 
gap) [24] to that of the covalently bonded diamond structure indicates a tetravalent behavior 
of Al atom. Thus, based on these observations, one can expect that the behavior of Al in 
Al n Li„ clusters would be similar to that of the tetravalent atom Si [25], Ge [26], Sn, and 
Pb [27]. Indeed, our present work shows that in certain clusters, structure and bonding of Al n 
is similar to that of the group IV A clusters. Moreover, earlier work on the heterogeneous 
Al-Li clusters has focussed on specific aspects like systematics in the geometry, stability, 
shell closure, magic behavior etc. Bonding in such clusters was discussed on the basis of 
total charge density as well as the difference of the self-consistent charge density, p sc f, and 
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the superimposed atomic charge densities p superimposed of the constituent atoms. On the 
contrary, analysis of the molecular orbitals (MOs) has revealed some unusual properties 
like aromaticity and antiaromaticity in all-metal Al-Li clusters. In our earlier work, we 
have examined some of the issues concerning the geometry and the stability of various Al- 
Li clusters viz. Al n Li n (n =, 1-10, 13) [28], Al n Li 7 (n=l-7) [29], AlLi n (n=l,8) [29] and 
Ali 3 Li„ (n= 1-4, 10, 19,20, 21) [30]. We have found that the geometries of these clusters were 
dictated by the geometry of the core Al cluster, enclosed in the Li cage. However, since 
these calculations were performed by the density based molecular dynamics (DBMD), an 
analyses of the bonding and electronic structure could not be done. In the present work, 
we study the systematics of the geometries, energetics, electronic structure and the bonding 
properties in various Al-Li clusters viz. Al n Li„ (n=l-ll), ALT, Al| — , Al 2 Li, Al 2 Li~, and 
Al 6 Lig using Kohn-Sham formulation of the density functional theory (DFT) within the 
pseudopotential and the generalized gradient approximation (GGA). The bonding in these 
clusters is analyzed via the electron localization function (ELF) [31] and the molecular 
orbitals. In section II we shall describe the computational details, followed by a discussion 
of the results in section III. 



II. COMPUTATIONAL DETAILS 

The ground state geometries as well as other low-lying structures were obtained in two 
stages. In the first stage, Born-Oppenheimer molecular dynamics (BOMD) [32] based on 
Kohn-Sham (KS) [33] formulation of density functional theory (DFT) was employed. The 
computer code used for this purpose was developed in our own group. The total en- 
ergy during each of the molecular dynamics step was minimized using damped equation 
of motion [34]. The calculations were performed using norm- conserving pseudopotentials of 
Bachelet [35] et al, in the Kleinman-Bylander [36] form with s part treated as nonlocal. The 
exchange-correlation potential was calculated using the local density approximation (LDA) 
given by Ceperley-Alder [37]. A cubic supercell of length 40a. u. with an energy cutoff of 
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~ 17 Rydberg was found to provide sufficient convergence of the total energy. Increasing 
the energy cutoff did not improve the total electronic energy significantly. 

The damped equation of motion scheme [34] permits use of a fairly moderate time step 
palOOa.u. Starting from a random configuration, the Al n Li n clusters were heated to 1400- 
1 700.fr, and allowed to span the phase space for a few thousand iterations. In our previous 
investigations [28-30] on Al n Li n clusters, it was found that the Li atoms segregate at the 
surface with the Al atoms forming an inner core. In order to avoid any such segregation 
during the molecular dynamics run, we have interchanged the Al and the Li atoms. This 
has ensured that the cluster visits all its local minima as well as the global minimum. At 
the end of each ionic displacement, the norm of the eigenstates defined as | hipi — | 2 
(where q is an eigenvalue corresponding to the eigenstate ipi of the Kohn-Sham Hamiltonian 
h) was within the range of 10~ 4 — 10~ 7 a.u. 

In the second stage, various low-lying structures were obtained by conjugate gradient 
and/or steepest descent [32] method starting from various suitable configurations during 
the molecular dynamics run. These configurations were then optimized using ultrasoft pseu- 
dopotentials [38] within the generalized gradient approximation (GGA) implemented in the 
VASP [39] package. The Perdew-Wang [40] exchange-correlation potential for GGA has 
been used. The size of the simulation cell was varied according to the cluster studied (see 
Table-I). 

Table - I 

Size of the supercell (in A) for various clusters, where n is the total number of atoms in the 

clusters. 
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n 


Simulation Cell (A) 


2,4,6 


16x16x16 


8,10,12 


18x18x18 


14,16.18,20 


20x20x20 



The geometries were optimized with a kinetic energy cutoff of the order of ~ 12 Rydbergs. 
This energy cutoff was sufficient for the convergence of the total energy. Increasing the 
energy cutoff did not improve the total electronic energy significantly. The structures were 
considered to have converged when the forces on each ion was less than O.OleV/A with a 
convergence in the total energy within the range of 10 -4 — 10 _6 eV. In general, we find that 
there are many isomeric structures nearly degenerate to the lowest energy state. Few of the 
structures can be obtained by interchanging the Al and the Li atoms or by rearranging the 
position of Li atoms. In the present work, we discuss only a few geometries corresponding 
to the lowest energies. 

The nature of the bonding has been investigated using the electron localization function 
(ELF) [31] along with the molecular orbitals (MO) as well as the difference charge density. 
The ELFs have been found to be useful for elucidating the bonding characteristics of a 
variety of systems, especially in conjunction with the charge density. The value of the ELF 
lies between and 1, where 1 represents a perfect localization of the valence charge density. 
The difference charge density is the difference between the self-consistent charge density, 
p sc f, and the superimposed atomic charge densities, p S uperim P osed, of the constituent atoms. 
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III. RESULTS AND DISCUSSION 



A. Structure, energetics and stability 

The ground state geometries and some of the low energy structures of Al n Li n , and 
Al 6 Li 8 clusters are shown in figures 1 (n—2-6, and Al 6 Li 8 ), and 2 (n=7-ll). The ground 
state geometries are divided into two broad categories, for (i) n—2-A (figures la(i)-lc(i)), 
and (ii) n=7-9, and 11 (figures 2a(i)-2c(i), 2e(i)), depending on the growth pattern. In the 
first category a quinted roof-like structure of Al2Li2 is seen (a quinted roof structure is a 
rhombus with a bend at its diagonal). The second category shows a pentagonal ring which 
is a signature of an icosahedral growth. The ground state geometries of n=5, 6, and 10 
clusters (figures Id, le, 2d) show neither a quinted roof nor a pentagonal ring. In what 
follows, we discuss the structure and stability of these clusters. The stability is discussed via 
the binding energy (E b ), the dissociation energy (AE) and the second difference in the total 
energy with respect to a single Al-Li pair (A 2 E), and the energy gap between the highest 
occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). 
These quantities are plotted in figure 3a (E 6 ), 3b (AE, A 2 E), and 3c (HOMO-LUMO gap) 
as a function of the number of Al atoms in the cluster and are defined as 

E 6 [Al n Li n ] = [ - E[Al n Li n ] + n (E[A1] + E[Li]) ] jln 
AE[Al n Li n ] = [E[Al n Li n ] - E[Al n _ 1 Li n _ 1 ]j - E[A1 Li] 
A 2 E[Al n Li n ] = [EfAL^LL^] - 2E[Al n Li n ]] - E[Al n+1 Li n+1 ] 

Let us recall that a pronounced maximum in the dissociation energy AE, along with a 
corresponding minimum in A 2 E, signifies stability. In addition, in figure 3d, we have plotted 
the smallest bond distances of Al-Al, Al-Li and Li-Li as a function of the number of Al atoms 
in the cluster. 

The ground state geometry of Al 2 Li 2 is a bent rhombus (figure la(i)). The planar 
geometry (figure la(ii)) is the low-lying structure of the Al 2 Li 2 . This cluster is one of the 
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most stable clusters of Al n Li n . A substantial rise in the binding energy (figure 3a) is seen 
when two Al-Li dimers combine to form Al 2 Li 2 cluster. Moreover, a peak in the dissociation 
energy (AE) (figure 3b), a minima in the second difference in energy (A 2 E) (figure 3b) and 
a rather large HOMO-LUMO gap (figure 3c) also signifies a relatively high stability of this 
cluster. Interestingly, the planar Al 2 Li 2 structure, with some bending, has been observed 
by Kumar [21] in the magic cluster Ali Li 8 . However, because of his restricted nature of 
search, he did not find this 8 valence electron system to have particular stability. 

A similar quinted roof structure is also seen in the ground state structures of the clusters 
Al 3 Li 3 and Al 4 Li 4 (figures lb(i) and lc(i)). However, the Al 3 Li 3 cluster is not energetically 
as stable as Al 2 Li 2 , as seen from figures 3a, 3b, and 3c. The lowest energy structure of 
Al 4 Li 4 shows two units of Al 2 Li 2 (figure lc(i)). These two units are arranged so as to form 
a rectangle out of the four Al atoms. Indeed, it is this rectangular shape that gives rise to 
an antiaromatic nature in this cluster [9,12]. This structure has been extensively studied by 
us [9,11,41] and others [12]. The two low-lying structures of Al 4 Li 4 (figure lc(ii) and lc(iii)) 
show a square-plane and a quinted roof of four Al atoms, respectively. 

In Al 5 Li 5 cluster, a square-plane of Al 4 is seen, with the fifth Al atom capping this plane 
forming a pyramid (figure Id). Such a square-planar geometry of Al 4 has been observed in 
the all-metal aromatic compounds Al 4 ~ and MAlJ (M = Li, Na, and Cu) [8]. Al 5 Li 5 shows 
the first three dimensional structure of Al n in Al n Li n clusters. Recall that the appearance of 
such three dimensional geometry in pure Al clusters [16] is seen for Al 6 . The Al 5 Li 5 cluster 
is the most stable among the clusters studied. This is clear from the binding energy, AE, 
A 2 E, and the HOMO-LUMO gap plots shown in figures 3a-3c, respectively. An interesting 
aspect of the Al 5 Li 5 cluster is the formation of a plane composed of four Li atoms with 
one Al atom at the center, akin to the face of a fee cell. Similar structure for Al 5 Nas was 
predicted by Dhavale [18] et al. A complete fee structure can be formed by fourteen atoms: 
six Al atoms at the face centered sites and eight Li atoms at the vertices of the cube. Such 
a structure was first studied by Shah [28] et al. using DBMD. They found this structure 
to be the lowest energy structure of Al 6 Li 8 . However, our investigation shows that this fee 
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structure, with slight distortion, is 1.1 leV higher than that the most stable structure. The 
lowest energy structure and the fee structure of this cluster is shown in figure If (i) and 
lf(ii), respectively. It can be noted that these two structures are nearly similar, the main 
difference being the elongation of one of the Al 4 square-planes in the fee structure into a 
rectangle. It is interesting to note that the series Al| _ , Li-AlJ, Al 5 Li 5 , and Al 6 Li 8 shows a 
three dimensional octahedral growth of Al cluster surrounded by Li atoms. This results in 
an effective derealization of the electron density (see section III B). 

The Al 6 Li 6 cluster also shows a geometry similar to that of Al 5 Li 5 , with distortions in the 
AI4 square-plane (figure le(i)). In the first low energy structure (figure le(ii)), a tetrahedra 
of Al 4 is seen, while the other low lying structure (figure le(ii)) shows an octahedra of Al 6 . 
The octahedral structure is the ground state geometry of pure Al 6 cluster [16]. 

The appearance of a pentagonal ring, which is a precursor to the icosahedral growth, 
is seen for n = 7 (AI13 is a distorted icosahedron in its ground state [42]). The ground 
state geometry of Al 7 Li 7 shows a pentagonal bipyramid which is capped by the Li atoms 
(figure 2a(i)). Appearance of a similar pentagonal ring in pure Al clusters (in neutral as well 
as singly charged form) is seen for Al 9 [16] as well as in tetravalent atom clusters viz. Sn, and 
Pb [27] for seven atom clusters. Incidentally, these clusters are isoelectronic to Al 7 Li 7 with 
28 valence electrons. The low energy structure of Al 7 Li 7 (figure 2a(ii)) as well as the ground 
state geometries of n — 8, 9, and 11 clusters (figure 2b(i), 2c(i) and 2e(i), respectively) 
also shows such pentagonal rings. However, the lowest energy structure of Al 10 Li 10 cluster 
(figure 2d(i)) is completely different from that of the Al 7 Li 7 to Al 9 Li 9 and AlnLin, in that, 
no signature of the pentagonal ring is seen. In fact, the geometry is similar to that of Sn 10 , 
i.e. a tetracapped trigonal prism (TTP) [43] with some distortions. This TTP structure 
is enclosed in a Li cage resulting in its distortion. It also shows magical behavior, since it 
exhibits a peak in the dissociation energy (AE), a minima in the second difference in energy 
(A 2 E), and a large HOMO-LUMO gap (figure 3a-3c, respectively). 

We have performed an analysis of the interatomic bond distances in order to understand 
the mixing and segregation behavior of the Al and the Li atoms. In figure 3d, we plot the 
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smallest bond distances of Al-Al, Al-Li and Li-Li as a function of number of Al atoms in 
the Al n Li„ cluster. We find that in smaller clusters (n < 4) the Al-Al bond is dominant. 
However, from n — 5 onwards, the Al-Li bond distance is comparable to that of Al-Al. 
Recall that at n — 5, the first three dimensional Al n structure emerges, thereby increasing 
the surface area. This helps the Li atoms to maximize the number of Al-Li bonds and 
thereby move closer to the Al atoms. 

It is well known that the density based method is not as accurate as the Kohn-Sham 
method, nonetheless, it has been extensively used to investigate the geometry, and stability 
of various Al-Li clusters [28-30] . Here, we make some pertinent comments on the difference 
in the geometry of the Al„Li n clusters obtained by these two methods. Our results show 
that the geometries of some of the clusters were in good agreement with those predicted 
by DBMD. For instance, the structure of Al 2 Li 2 was found to be a quinted roof by both 
the methods (see ref. 28 for DBMD results). The geometries of other clusters like Al 7 Li 7 to 
Al 9 Li 9 show some substitutional disorder. However, there are certain clusters where the 
geometries were quite different by the KS method. For instance, the structures of, Al 3 Li 3 , 
predicted by DBMD, is a capped trigonal bipyramid [28]. When optimized by KS, this 
structure gets Jahn- Teller distorted into a structure with a quinted roof of Al 2 Li 2 . The 
evolution in the geometry by DBMD method shows an earlier appearance of three dimen- 
sional structure of Al n at n — 4 (Al 4 is tetrahedra in Al 4 Li 4 ), as compared to the pyramidal 
structure of Al 5 in Al 5 Li 5 by KS method. 

B. Bonding 

1. Al 2 Li 2 

As noted above, Al 2 Li 2 is one of the most stable cluster. This cluster has 8 valence 
electrons, which corresponds to a closed shell in the jellium sense. However, the behavior 
of the MOs as well as the eigenvalue spectrum does not resemble to that of the SJM. 
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Kumar [21] found that the bonding between Al-Li in the planar Al 2 Li 2 structure is ionic. 
On the contrary, our investigation on the bent rhombus structure does not show signature 
of any significant charge transfer from the Li atoms to the Al atoms. However, an analysis 
of the difference charge density shows that, a charge transfer from the Li atoms to the 
bonding region between the Al and Li atoms takes place (figure not shown). In order to get 
better insight of the bonding in Al 2 Li 2 , we have studied various clusters viz. Al 2 , Al 2 , Al 2 ~, 
Al 2 Li, Al 2 Li~, and Al 2 Li 2 . This series represent a growth of Al 2 Li 2 from the Al 2 dimer. In 
figure 4, we show the eigenvalue spectrum, and the Al-Al bond distance for the optimized 
geometries of these clusters. The eigenvalue spectrum of Al 2 (figure 4a), in neutral, singly 
and doubly charged form, shows a triply degenerate HOMO state. These states are partially 
filled two 7r and one a bond (figure not shown). Successive substitution of each electron in 
the charged dimer, by Li atoms splits the HOMO state reducing the degeneracy. The first 
Li atom forms a a bond with the p y — p y 7r orbital of the Al 2 dimer (isodensity surfaces not 
shown). The second Li atom further splits the HOMO state introducing a substantial gap of 
0.68eV. However, in this case, a hybridization of the s orbital of Li atom and the 7r orbital 
of Al 2 takes place. In figure 5, we show the second HOMO and the HOMO state of the 
Al 2 Li 2 cluster. Clearly, these states shows the 7r bonded orbitals of Al 2 hybridized with the 
s orbitals of the Li atoms. In figure 4b, we plot the Al-Al bond distance in these clusters. 
A contraction of the Al-Al bond upon addition of an electron or a Li atom is seen. The 
contraction of the bond on addition of a Li atom being more than that due to the addition 
of an electron. Thus, the hybridization of the s orbitals of Li with the n orbitals of Al 2 
enhances its stability. 



2. AkLi A 

The bonding in Al 4 Li 4 cluster has been extensively studied by us [9,11,41] and oth- 
ers [12]. This cluster has been found to be antiaromatic with four 7r-electrons. For the 
sake of completeness, we discuss the bonding in this cluster. This cluster, as discussed ear- 
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lier, is composed of two Al 2 Li 2 units arranged edge-to-edge to form a Al 4 rectangle. This 
rectangular structure leads to antiaromaticity. In this cluster sp 2 hybridization of Al takes 
place leaving one empty unhybridized p orbital. The valence electron of each of the four Li 
atoms is then transfered to this empty p z orbital, thus providing four 7r-electrons for antiaro- 
maticity. Interestingly, its first low-lying structure (figure lc(ii)) can also be considered as a 
candidate for antiaromaticity. An analysis of the MOs shows that this cluster also has four 
7r-electrons similar to that of the lowest energy structure, thereby conforming the structural 
and electron count criteria for antiaromaticity. However, a detailed analysis of the magnetic 
field induced ring currents has to be done in order to understand the antiaromaticity in this 
structure. 



3. Al^L% and Al^Lig 

As discussed in section III A, the structures of Al 5 Li 5 and Al 6 Li 8 (in fee geometry) 
clusters show a square-plane of Al 4 . Such Al 4 square-plane in Al|~ and Li-AlJ has lead to 
aromaticity in these clusters [8]. Moreover, these clusters (All - , Li-AlJ, Al 5 Li 5 , and Al 6 Li 8 ) 
show a three dimensional octahedral growth of the Al cluster. It is, therefore, interesting to 
discuss the change in bonding, and hence the effect on aromaticity in these clusters. Recall 
that Al| _ and Li-AlJ are aromatic clusters with two 7r-electrons [8]. The HOMO state in 
these clusters are completely delocalized 7r orbitals. There are two more delocalized a bonds: 
one composed of radial p orbitals, and the other composed of perpendicular p orbitals. Due 
to the presence of a Al 4 square-plane, the bonding in Al 5 Li 5 and Al 6 Li 8 clusters is expected 
to be similar to that in Al|~. In that case, these clusters would have two 7r-electrons similar 
to that of Al|~ and Li-AlJ, and show some signature of aromaticity. In figure 7a and 7b, we 
show the HOMO states of Al 5 Li 5 and Al 6 Li 8 . Interestingly, the HOMO state of Al 5 Li 5 does 
show such a delocalized n bond between the Al 4 square (figure 7a). However, a lone electron 
is seen at the fifth Al atom that caps this square. This reduces the number of n electrons 
count to one, thereby losing the aromaticity. The HOMO state of Al 6 Li 8 , on the other hand, 
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is quite different than that seen in the previous clusters (see figure 7b). It shows a localized 
bond composed of p orbitals of the Al atoms. In figure 6, we show the eigenvalue spectrum 
of Al 6 Li 8 , along with that of the other Al n Li„ (n =1-11) and Sni clusters. It is clear from 
this that the eigenvalue spectrum of Al 6 Li 8 shows nearly a jellium-like structure. However, 
an analysis of the behavior of the molecular orbitals (figure not shown) of Al 6 Li 8 shows that 
the lower six states are jellium-like, whereas, the higher states show localized bonds arising 
out of p orbitals of the Al atoms. The isodensity surfaces for the 6 th state of Al 5 Li 5 and the 
10 th state of Al 6 Li 8 are shown in figure 7c, and 7d, respectively. These are a-bonded states, 
composed of radial p orbitals and are analogous to the radial p orbitals of Al 4 ~, as noted 
by Wang et al [8]. Since, these bonds involve all the Al atoms, they exhibit a substantial 
derealization. 



4. A^Lin (n=7-ll) 

The bonding in these clusters (Al n Li n , n = 7—11), is quite different from that of 
the smaller clusters. The analysis of the behavior of the MOs shows that the eigenvalue 
spectrum can be classified into two groups: the lower n states are jellium-like, whereas, the 
higher states forms a complex band of localized bonds arising out of the interaction of the p 
electrons of each of the Al atoms. The lower states can be identified as Is, Ip, Id, 2s,..., in 
conformation to the SJM. Typical representative isodensity plots of these states are shown 
in figure 8. They represent the Is, Ip, Id and 2s states of the jellium composed of the Al-3s 
orbitals. In figure 6, we show the eigenvalue spectrum of Al n Li n (n =1-11) and Sni clusters. 
This figure shows a clear energy separation of lower jellium-like n states of the clusters with 
n > 6, from the higher orbitals. In order to study this more extensively, we have performed 
a spherical harmonics analysis of the KS orbitals [44]. We find that the lower n orbitals 
have a predominant s character (75-85%) indicating that these state are arising out of the 
3s orbitals of the Al atoms, which hybridizes to form a jellium. However, this description 
cannot be extended to the higher states. The higher states are composed of A\-3p orbitals. 
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The localization characteristics in the bonding in the Al-Li clusters is also analyzed via 
the electron localization function plots. Previous investigations show that even in large 
Al-Li clusters, the Al-Al bond is covalent. However, our earlier work on the structural 
properties shows a clustering of Al atoms enclosed in Li cage. This clustering could lead to 
a derealization of the electron density. In order to study this, we have analyzed the electron 
localization function in the Al n Li n clusters. In figure 9, we show the isovalued surface of 
the ELF for the clusters Al 2 Li 2 , Al 5 Li 5 , and Ali Lii , at the values 0.85, 0.8, and 0.72, 
respectively. These plots shows a localization of the electron density along the Al-Al bond 
(increasing the value of ELFs does not show any bond between Al-Al). It has been noted 
by Silvi et al [31], that an ELF value of about 0.7 or greater is a indication of a localized 
bond in that region. Thus, the Al-Al bonds in these clusters are predominantly covalent, in 
accordance with the previous observation. However, as the cluster grows in size, the degree 
of localization of the electron density in these clusters is reduced. Further, the localization 
in Al 5 Li 5 is only in the plane forming the square (for out of plane, the Al-Al bond is seen 
for the ELF value of 0.7). A similar behavior of the ELF for Al 7 Li 7 (not shown in figure) is 
seen, in that, a localized bond along the pentagon (ELF value = 0.78) and delocalized bond 
in the perpendicular direction (ELF value = 0.6) is seen. 

IV. TETRAVALENT BEHAVIOR OF Al IN Al n Li n CLUSTERS 

Previous investigations on mixed Al-Li clusters [21,45,46], have shown a substantial 
charge transfer from the Li atom to the Al atom. Due to such charge transfer, Al be- 
haves as a tetravalent atom in mixed Al-Li clusters with nearly 50:50 percent concentration. 
In this section, we present the evidence of such behavior. The analysis of the geometries 
of Al n in Al n Li„ shows remarkable similarities to those of the clusters of the tetravalent 
atoms viz., Si, Ge, Sn and Pb. The transition from planar to three dimensional geometry 
at n — 5, the appearance of a pentagonal ring at n = 7, and the formation of the distorted 
TTP structure of 10 atom cluster are some of the features of the clusters of tetravalent 
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atoms. Such behavior of Al„ in Al n Li„ clusters is also observed. Moreover, the eigenvalue 
spectrum, in figure 6, of the Al 10 Li 10 show remarkable similarity to that of the Sn 10 cluster. 
The bonding between Al-Al in these clusters is covalent, which is similar to that in the 
clusters of tetravalent atoms. Finally, the resemblance of the DOS of the Al-Li in its most 
stable bulk phase i.e. B32, to that of diamond [24] also indicates a tetravalent behavior of 
the Al atom in Al-Li systems with 50:50 concentration. 

V. CONCLUSION 

In the present work, we have reported the systematic investigation of the geometries, 
energetics, electronic structure and bonding in various Al-Li clusters viz. Al n Li n (n=l- 
11), AL~, Al 2 ~, Al 2 Li, Al 2 Li~, and Al 6 Li 8 using BOMD within the framework of density 
functional theory. The bonding in these clusters was discussed via the electron localization 
function as well as the behavior of the molecular orbitals. The growth pattern is divided 
into two broad categories: first consisting of a bent rhombus of Al 2 Li 2 (n=2-4) and second 
consisting of a pentagonal ring (n=7-9,ll). We find that the 8, 20, and 40 valence electron 
systems to be magic, exhibiting a peak in the dissociation energy, a minima in the second 
difference in energy, and a large HOMO-LUMO gap. The structural transition of Al n in 
Al n Li„, from two dimension to three dimension increases the surface area, thereby helping 
the Li atoms to maximize the number of Al-Li bonds, as a result of which the Li atoms 
move closer to the Al atoms. In clusters with more than six Al atoms, a charge transfer 
from Li atoms to Al atoms makes Al behave as a tetravalent atom like Si, Ge, Sn, and Pb in 
Al n Li„ clusters. These negatively charged Al„ structure is then stabilized by the positive Li 
environment. The behavior of the MOs in these clusters (n > 6) can be divided two groups: 
a lower group of jellium-like states arising out of the s electrons of the Al atoms and a higher 
group of localized bonds formed through the interaction of the p electrons on each of the 
Al atoms. The formation of the three dimensional Al n structure destroys the aromatic and 
antiaromatic nature in these clusters. 
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List of Figures 

1 The ground state and the low-lying geometries of the Al„Li n (n—2-6) and Al 6 Li 8 clus- 
ters. The black circles represent Al atoms and the white circles represent the Li atoms. 
The lowest energy structure is represented by (i). 

2 The ground state and the low- lying geometries of the Al n Li„ (n=7-ll) clusters. The 
black circle represents Al atoms and the white circles represent the Li atoms. The 
lowest energy structure is represented by (i). 

3 (a) The binding energy per atom (in eV/atom) of the Al n Li„ cluster (n=l-ll). (b) The 
dissociation energy (AE) and the second difference in the total energy (A 2 E) of the 
Al n Li„ clusters (n=2-10) with respect to a Al-Li dimer. (c) The HOMO-LUMO gap 
(in eV) of the Al n Li„ cluster (n=l-ll). (d) The minimum interatomic distance Al-Al, 
Al-Li and Li-Li (in A) of the Al n Li n cluster (ra=l-ll). 

4 (a) The eigenvalue spectrum (in eV). The continuous lines correspond to occupied 
states and dashed lines correspond to the empty states. The numbers on the right 
indicate the degeneracy of the states. The numbers on the left indicate the occupancies 
of the states. All other occupied states have two electrons, (b) The Al-Al bond distance 
(in A). 

5 The isodensity surface of the second HOMO and the HOMO state of Al 2 Li 2 at l/5 th of 
its maximum value. The black circles represent Al atoms and the gray circles represent 
the Li atoms. 

6 The eigenvalue spectrum (in eV) of Al„Li n (n=l-ll), Al 6 Li 8 , and Sn 10 clusters. The 
continuous lines correspond to the occupied states and the dotted lines correspond to 
the empty states. The numbers on right indicate the degeneracy of the states. All 
states are doubly occupied. The thick dashed line shows the n th state of these clusters, 
where n is the number of Al atoms. For Snio the thick dashed line corresponds to the 
10 th state. 
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7 (a) The isodensity surface of the HOMO state of Al 5 Li 5 , at l/5 th of its maximum. 

(b) The isodensity surface of the HOMO state of Al 6 Li 8 , at l/5 th of its maximum. 

(c) The isodensity surface of the 6 th state of Al 5 Li 5 , at l/5 th of its maximum, (d) The 
isodensity surface of the 10 th state of Al 6 Li 8 , at l/8 th of its maximum. In all the figures, 
the black circles represent Al atoms and the gray circles represent the Li atoms. 

8 The isodensity surface of various orbitals. These are: (a) 1 st orbital of Al 7 Li 7 at l/5 th 
of its maximum value, (b) 3 rd orbital of AlgLig at l/5 th of its maximum value, (c) 8 th 
orbital of Al 9 Li 9 at l/10 th of its maximum value, and (d) 10 th orbital of Ali Lii at 
l/5 th of its maximum value, representing the Is, lp, Id, and 2s state of the jellium 
composed of the Al-3s orbital. In all the figures, the black circles represent Al atoms 
and the gray circles represent the Li atoms. 

9 The isodensity surface of the electron localization function (ELF) of: (a) Al 2 Li 2 at 
the value 0.85, (b) Al 5 Li 5 at the value 0.80, and (c) Al 10 Li 10 at the value 0.72. In all 
the figures, the black circles represent Al atoms and the gray circles represent the Li 
atoms. 
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